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We uncover the general mechanism producing the dark energy. This is only based on well known 
quantum physics and cosmology. We show that the observed dark energy originates from the 
cosmological quantum vacuum of light particles which provides a continuous energy distribution 
able to reproduce the data. Bosons give positive contributions to the dark energy while fermions 
yield negative contributions. As usual in field theory, ultraviolet divergences are subtracted from 
the physical quantities. The subtractions respect the symmetries of the theory and we normalize 
the physical quantities to be zero for the Minkowski vacuum. The resulting finite contributions to 
the energy density and the pressure from the quantum vacuum grow as log a(t) where a(t) is the 
scale factor, while the particle contributions dilute as l/a 3 (t), as it must be for massive particles. 
The dark energy equation of state P — w(z) Tt turns to be w(z) < — 1 with w(z) asymptotically 
reaching the value —1 from below. A scalar particle can produce the observed dark energy through 
its quantum cosmological vacuum provided: (i) its mass is of the order of f 0~ 3 eV = 1 meV, (ii) 
it is very weakly coupled and (iii) it is stable on the time scale of the age of the universe. The 
axion vacuum thus appears as a natural candidate. The neutrino vacuum (especially the lightest 
mass eigenstate) can give negative contributions to the dark energy. We find that w(z — 0) is 
slightly below —1 by an amount ranging from — f .5 x 10~ 3 to —8 x 10~ 3 while the axion mass results 
between 4 and 5 meV. We find that the universe will expand in the future faster than the de Sitter 
universe, as an exponential in the square of the cosmic time. Dark energy arises from the quantum 
vacua of light particles in FRW cosmological space-time in an analogous way to the Casimir effect 
in Minkowski space-time with non-trivial boundaries. 
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I. INTRODUCTION AND RESULTS 



Since the discovery of the dark energy in the present universe [l| an intense observational activity has improved 
our knowledge about it @, [|[ . Many theoretical models have been proposed to explain the dark energy. A common 
feature of such proposals is that all of them postulate extra physics as source of the dark energy 0, [j| . 

We show in this article that the dark energy present today in the universe originates from the cosmological quantum 
vacuum of light particles. This is a vacuum effect which unavoidablely appears when quantum fields evolve in a 
cosmological space-time. That is, dark energy is generated by a mechanism based on well known quantum physics 
and cosmology. Bosons yield positive contributions to the dark energy while fermions give negative contributions. 

We find that the scale of the contributions to the dark energy is of the order of 

M 4 

2(4^ l0gZdcc ' 

where M is the particle mass and z& ec is the redshift when it decoupled from the early universe plasma. 

Generally speaking, the energy of a quantum field is the sum of the vacuum contribution plus particle contributions. 
It is known that the vacuum energy of a quantum field dissipates into particles when the field evolves coupled to other 
fields or to itself [I(|-[l3j]. Dissipation into fermions is reduced by Pauli blocking [ll|, [l2j]. Electrons, protons and 
photons are coupled to photons and therefore, their vacuum energy dissipates through photon production well before 
recombination, that is, when the temperature of the universe was 1 MeV or more. Unstable particles cannot produce 
long-lasting vacuum effects. Only a very weakly coupled stable particle can produce a vacuum energy contribution 
lasting for times of the order of the age of the universe, that is, a vacuum energy contribution measurable today. 

Since the dark energy is known to be positive, bosons must dominate the cosmological vacuum energy. The scale 
of the boson mass must be in the meV range since the observed dark energy density has the value [25(, 

/?A = Pc = (2.39 meV) 4 , 1 meV = 1CT 3 eV . (1.1) 

Spontaneous symmetry breaking of continuous symmetries is a natural way to produce massless scalars (Goldstone 
bosons) in particle physics. Furthermore, a slight violation of the corresponding symmetr y can g ive a small mass to 
such scalar particle. Axions, majorons and familons have been proposed in these grounds [161 1 X 71. |21|. 
In addition, the lightest neutrino can give a negative contribution to the dark energy. 

Neutrinos in the universe are known to be free for temperatures T < 1 Mev which corresponds to redshifts z < 
6 x 10 9 @, [HI]. That is, we can describe their evolution as free fermions in the cosmological FRW universe. 

Axions with masses M ~ 1 meV are free for temperatures T < 10 6 GeV which corresponds to redshifts z < 10 19 18]. 
They can be considered as free scalars in the cosmological FRW universe. Both, the axion and neutrino decoupling 
happens during the radiation dominated era. Before decoupling, the non-negligeable interaction of the corresponding 
particles made dissipation important and therefore the vacuum energy can only become significant after decoupling. 
Therefore, we can restrict ourselves to study the free quantum field evolution in the cosmological space-time after 
decoupling. 

We investigate the evolution of scalars and fermions as an initial value problem (Cauchy problem) for the corre- 
sponding quantum fields on a cosmological space-time. We find that the initial temperature has a negligible effect on 
the vacuum energy for late times. 

Both axions and neutrinos can lead to vacuum effects lasting cosmological time scales. Any of the two heavier 
neutrino mass eigenstates v-i and v% would produce a large negative dark energy in the (50 meV) 4 range. Hence, 
either the heavier neutrinos V2 and annihilate with their respective anti-neutrinos in a time scale of the age of 
the universe or a stable scalar particle with mass in the > 50 meV range must be present in order to reproduce 
the observed value of the dark energy eq. p.ip . However, we find that the second possibility is inconsistent with the 
observed dark energy equation of state. 

An effective four fermions interaction with strength characterized by M'~ 2 where M' is a mass scale can make the 
heavier neutrinos unstable. The mass scale M' should be M' < 1 MeV or M' < 10 MeV for the direct and inverse 
neutrino mass hierarchies. 
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In conformal time, the scalar and fcrmion fields rescaled by the scale factor turn out to obey equations of motion 
similar to Minkowski space-time with time-dependent masses 



M 2 a 2 ( V ) - 



a"(v) 

0,(1]) 



X(x,rj) = , 



[i ft — to a(rj)] ip(x, rj) = 



(1.2) 



Here, x and ip are, respectively, rescaled scalar and fermion fields, V 2 is the usual flat space laplacian and iff) is the 
usual Dirac differential operator in Minkowski space-time in terms of flat space-time Dirac matrices. 

There are two widely separate scales in the field evolution in the cosmological space-time. The fast scale is the 
microscopic quantum evolution scale, typically ~ 1/M ~ 1/tti, where M and to are the scalar and fermion masses, 
respectively. The slow scale is the Hubble scale \/H of the universe expansion. When M ~ m ^> H, M 2 ^> 
a" (rf) I a 3 (77) , and hence the scale factor can be considered as constant. Therefore, for the fields x and V 7 the cosmological 
quantum field evolution is just the Minkowski evolution with effective masses M 2 a 2 and to a, respectively, as seen 
from eg- CCH) . 

Energy density, pressure, and field density express in field theory as products of field operators and their derivatives 
at equal space-time points. Such expressions are ultraviolet divergent and need to be subtracted. The subtractions 
respect the symmetries of the theory and we normalize them such that the physical quantities are zero for the vacuum 
in Minkowski space-time. The finite resulting quantities grow as log 0,(77). This is analogous to the high-energy growth 
of renormalized one-loop Feynman graphs. 

That is, the energy density and the pressure get contributions from the quantum vacuum that grow as log 0.(77) 
while the particle contributions dilute as l/a 3 (77), as it must be for massive particles. 

We obtain for the vacuum energy density and pressure of scalar and fermion fields with mass M and to, respectively, 



< n > (77) 



< p > (v) 



a(7j)»a d cs,<Mcf 



M 4 



a(!j)»a d cs,adcf 



2(4tt) 2 

M 4 



log 0(77) +65 



1 



77i 



(4tt) 5 



log 0(77) + hi 



2(4tt) 



log 0(77) + b s + — 



(4tt)< 



log 0.(77) +°F + ^ 



(1.3) 
(1.4) 



Where b$ and bp take into account the initial values of the scale factor Od cs and Odcf (at the decoupling time) of the 
scalars and fermions, respectively. H = 1 for Majorana fermions and H = 2 for Dirac fermions. 



We therefore obtain for the equation of state, 

< P > (V) a.(r,)»a dca ,a dc 



w(r)) 



<H>( V ) 



2 N 



(1.5) 



log 0(77) - 



2 Km 4 

TP 



That is, we find w(rj) < —1 with 10(77) asymptotically reaching the value —1 from below. 

It is convenient to express the scale factor in terms of the redshift. Taking into account that bs and bp contain the 
initial values of the scale factor yields, 



0(77) e b 



l + z s 
1 + z 



a(v) e b 



l + Zp 



where z$ (zf) is the redshift when the scalar (fermion) field decoupled. For neutrinos we know that zf 
while z$ ~ 2.2 x 10 18 for axions with mass ~ 1 meV |18| . 
We find from eqs.JOJ) and fTB]) . 



< H > (z) 



1 



< P > (z) = - 



2 (An) 2 
1 



2(4tt) 



M 4 logz s - 2 H to 4 log z F - (M 4 - 2 H to 4 ) 
M 4 \ogz s - 2 H to 4 logz_F - (M 4 - 2 N to 4 ) 



log(l + z) + 



log(l + z) - — 



(1.6) 
6 x 10 9 @ 

(1.7) 



where we used that zs ^> 1, Zf ^> 1- 
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We identify the vacuum energy density today < TL > (z — 0) with the observed dark energy pa- We can then write 
eas.lfO]), (fT3|) and (UTTJ) as, 



PA 



2(4?r) 



M log z s 



2 K m ( log z F 



<H> (r)) = p A 



w(r]) + 1 

where ao is the scale factor today and 



a(?))>a dcB ,a dcf 



i + ^iog^M 

M 4 - 2 H m 4 



6(4^) 2 PA 1 + & log^g) 



/8m = 



2 H to 4 

M 4 



(1.8) 



(1.9) 



log 2 S 



1 2 H to 



4 r 



M 4 



log z F 



(1.10) 



That is, the vacuum energy density at late times after decoupling grows as the logarithm of the scale factor and the 
equation of state asymptotically approaches —1 from below. 



2 Km 4 , / 2ttm 4 

lo & z S T7I- lo S z F ~ 1 - 



The equation of state as a function of z takes the form 

1 

w(z) + 1 = 
For z = it becomes today: 

w(0) + 1 = 



1 - 



2 N to 4 
M 4 



M 4 



M 4 



2Hm 4 
M 4 



logz s 



1 2 H m 



M 4 



log z F 



log(l + z) 



M 4 - 2 H m 4 
6(4 7t) 2 pa ' 



(1.11) 



(1.12) 



The scalar and fermion masses are constrained by the value of the dark energy today eq. (jl.ip . This gives the positivity 
requirement: 

M > (2 N)3 m , 
as well as the expression for the mass of the scalar particle: 

10.1 mcV 



M = 



(l0gZ5-i)* 



.3.90 meV/ 



(1.13) 



The neutrino contribution to the dark energy can be ignored when to <C 1 meV and when the vacuum neutrino 
contribution dissipates in the time scale of the age of the universe as mentioned before. The mass of the lightest 
neutrino is not yet known (only neutrino mass differences are known). We will consider that the lightest neutrino 
mass is either m — 3.2 meV [2(| [27| or zero [2{| . 

More specifically, we set zs ~ 2.2 x 10 18 assuming the scalar to be an axion with mass ~ 1 meV in eqs. (|1.12j> and 
(|1.13[) and we get for the axion mass and for the equation of state today, 



3.96 meV < M < 4.66 meV , 
-0.00794 < w(0) + 1 < -0.00156 



(1.14) 



The left and right ends of the intervals in eq. (|1.14p correspond to no neutrino contribution and to the lightest neutrino 
as a Dirac fermion with mass m = 3.2 meV, respectively. 
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We see that w(0) is slightly below —1 by an amount ranging from —1.5 x 10~ 3 to —8 x 10~ 3 while the axion mass 
results between 4 and 5 meV which is within the range of axion masses allowed by astrophysical and cosmological 
constraints 

If the scalar particle is not the axion, the value of 3> 1 will depend on the dynamics of such scalar particle. 

In general, we express the contribution of the quantum vacuum of light particles to the dark energy and pressure 
in terms of two parameters: the particle masses and the redshifts when they decoupled. There is also a dependence 
on the number of states per particle (1 for a scalar, 2K for a fermion). 

We uncover in this paper the general mechanism producing the dark energy. This mechanism is only based on well 
known quantum physics and cosmology. The dark energy appears as a quantum vacuum effect only due to the 
cosmological space-time expansion. In Minkowski space-time no dark energy is induced. The dark energy arises for 
a quantum field in the cosmological context in an analogous way to the Casimir effect in Minkowski space-time with 
non-trivial boundary conditions in space. 

All physical (finite) results are independent of any energy cutoff as well as of the regularization method used. 

We obtain and solve the self-consistent Einstcin-Friedmann equation for the scale factor when the dark energy 
dominates and the universe expansion accelerates. The growth of the energy density eq. (|1.3|) as the logarithm of the 
scale factor implies an expansion faster than in de Sitter space-time. More precisely, we find that the universe will 
reach in the future an asymptotic phase where it expands as 



a(t) a(today) e ClH ° t + C2 



(1.15) 



where 



ci 



fl A = 0.87 



0.00452 < c 2 < 0.00872 



(1.16) 



and Ho stands for the Hubble parameter today. The left and right ends of the interval for C2 in eq. (|1.16|) correspond 
to no neutrino contribution and to the lightest neutrino being a Dirac fermion with mass m — 3.2 meV, respectively. 

Notice that the time scale of the accelerated expansion is huge ~ l/-ffo = 13-4 Gyr. The quadratic term dominates 
over the linear term in the exponent of eq. (|1.15|) by a time t ~ 100/ifo to 200/iJo- 

In this accelerated universe, we see from the Friedman equation and eq. t|l .3[) that the Hubble radius 1 /H decreases 
with time as 1/[Hq y/log a(t)]. 

This paper is organized as follows: in sec. II and III we review the dynamics of scalar and fermion fields on 
cosmological space-times, respectively. In sec. IV we find the main physical quantities and the equation of state for 
late times. The dark energy is found as a result of the cosmological quantum vacuum contributions from light particles 
in sec. V. We obtain the future self-consistent evolution of the universe in sec. VI. We discuss relevant related issues 
in sec. VII and present our conclusions in sec. VIII. An appendix is devoted to the equivalence between different 
regularization methods. 



II. SCALAR FIELDS IN COSMOLOGICAL SPACE-TIMES 



We consider a massive neutral scalar field tp in a FRW geometry defined by the invariant distance 

ds 2 = dt 2 - a 2 (t) dx 2 . 

The Lagrangian density is taken to be 



(2.1) 



2 v y 



It is convenient to use the conformal time r\ 



and the conformally rescaled field x(af, rj), 



dt 



(2.2) 



a(t) 



X(x,r)) = a(t) ip(x,t) . 



(2.3) 
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The action (after discarding surface terms that do not affect the equations of motion) reads: 



A 



M - \ I 



d 3 x dr\ 



X ,2 -(Vx) 2 -M 2 W x 2 



where primes denote derivatives with respect to the conformal time 77 and where 

a" (77) 



M. (rj) = M a 2 {r)) 



0(77) 



(2.4) 



(2.5) 



plays the role of an effective mass squared. Therefore, the rescalcd field x(x, 77) obeys the equation of motion, 

X" -V 2 x + M 2 (r ] )x = 0. (2.6) 

The evolution of x{x,r]) IS like that of a scalar field in Minkowski space-time with a time-dependent mass squared 
M 2 ^). 

The solution for the field y (x, t) can be Fourier expanded as follows, 



<p(x,rj) = 



1 



d 3 k 



0(77)7 (2^)3 2^0 



HMv) e tk - x + al <j>%(r,) e 



(2.7) 



where 



E = Jk^+M 2 , 



and Mi is the effective mass M.{rj) at the decoupling time (initial time) for the scalar field evolution. The mode 
functions <j>k(i]) obey the evolution equations, 



— + k 2 + M 2 a 2 (r])- 



a"(r>) 



a(rf) 



(2.8) 



We choose the initial state as the vacuum state which is here (at decoupling) a thermal equilibrium state at temperature 
T. However, as we see below [eq. (|4.10|) ]. the effect of the initial temperature on the vacuum energy is negligible for 
late times. The Fock vacuum state |0 > is annihilated by the operators aj. Therefore, we have as initial conditions 
for the mode functions, 



fo(0) = l , <f/ k (0) = -iEo 



(2.9) 



These initial conditions describe the Bunch-Davies vacuum when they are applied at asymptotically earlier times in 
the past (77 — > —00) 0,Q. See the discussion in sec. IVHI below. 

The time-dependent creation and annihilation operators obey the canonical commutation rules, 



2E (2tt) 3 5(k-k') 



(2.10) 



The energy- momentum tensor for a scalar field is given by 0] , 

Tpv = d^ip d v ip - - [d x ip d x tp - M 2 ip 2 ] . 

Its expectation value has the fluid form 

< T° > = < Us > (77) , <Ti> = -5\ <P S > (77) , 

since we consider homogeneous and isotropic quantum states and density matrices. In conformal time the hamiltonian 
density and the pressure take the form 



Hs(r)) 



j7j~T {[X'{x,ri) - a(7?) H( V ) x&v)? + (V X (x, i])? + a 2 (v) M 2 X 2 (x,v)} > 



2 a 4 (77) 



X'(x, 77) - 0(77) H(r\) x(x, v)} 2 + g (Vx) 2 (f, 77) 



(2.11) 



where H (77) stands for the Hubble parameter 

„, % d In aft) 1 da ,„ „. 

It is convenient to consider the conformal energy and pressure, 

s s (v) = a\ri) <H S > fa) , p S (v) = a A {r,) < P s > (r?) . (2.13) 

We find the trace of the energy-momentum tensor from eqs. (|2.11[) . 

a\v) ln s (v) - 3 P S (V)] = a 2 (v) M 2 X 2 - [(*' ~ a h X f ~ (V X ) 2 - a 2 ( V ) M 2 X 2 ] . (2.14) 

Ignoring the bracket term in the right hand side yields the virial theorem. Although this bracket term is nonzero, its 
space and time average is zero: 

' dr, [ d 3 x [( X ' -ah X ) 2 - (V X ) 2 - a 2 ( V ) M 2 X 2 } A>> ^ /M . 



A 



In addition, this bracket can be neglected for late times as we shall see below. 
Therefore, we have for the expectation values, 

£s(v) ~ 3 Ps (v) = M 2 a 2 ( v ) E s fa) - a 4 (r?) V( V ) (2.15) 

where 

Esfa) =< X 2 (x, V) >= a 2 (v) < f 2 (x, V) > , (2.16) 
and V stands for the expectation value of the virial 

V( V ) =<( X '-ah X ) 2 - (V X ) 2 - a 2 ( V ) M 2 X 2 > . 

Using the equations of motion (|2.6p we obtain for the time derivative of the energy density eq. (|2. 13|) , 

w = 5 m2 Ss(r/) _ a{v) m v{71) ■ (2 - 17) 

This relation in conformal time implies the usual continuity equation in cosmic time 

^<H S > +3H( V ) (< H s > + < Ps >) = . (2.18) 
at 

Thus, we see from eas. (|2.15[) and (|2.17[) that there is only one independent quantity among es(r)), Ps(v) an d ^s(v)- 

III. FERMION FIELDS IN COSMOLOGICAL SPACE-TIMES 

The Lagrangian density for fermions is taken to beQ] 

C = y/^^U^ -ml* . (3.1) 
The 7^ are the curved space-time Dirac 7 matrices and the fermionic covariant derivative is given by 

Z>„ = + \ h c ,7 d }V: (DM 
D^V*, = d^Vdu - r*, V dX 

where the vierbcin field is defined as 

g* v = v»v?rf b , 
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r] a b is the Minkowski space-time metric and the curved space-time matrices 7 M are given in terms of the Minkowski 
space-time ones 7° by (Greek indices refer to curved space time coordinates and Latin indices to the local Minkowski 
space time coordinates) 

7" = 7 a K M , {^,l v }=2g^. 

In conformal time the vierbeins V£ are particularly simple 

K P = a( V ) Si , (3.2) 

where a(rj) = a(t(rj)) is the scale factor as a function of the conformal time and we call a(r) = 0) = Ode- The Dirac 
Lagrangian density thus simplifies to the following expression 



-g * (i 7 M - m) * = a** [i ft — m 0(77)] (a 2 \E< 



(3.3) 



where if) is the usual Dirac differential operator in Minkowski space-time in terms of flat space time 7° matrices. 
Therefore, the Dirac equation in the FRW geometry is given by 



[i ft — m a(i~i)] a 2 ^{x,ri) 
The solution ^(x, 77) can be expanded in spinor mode functions as 







E 



d 3 k 



b^Ux&ri + dl^Vxi-k,!!) 



where 



e = \Jk 2 + m 2 a\ c , 

and the spinor mode functions U, V obey the Dirac equations 

i 7 d v — 7 • k — m a(r]) U\(k,rf) = 
i 7 drj + 7 • k — m a(rj) V\(k,rj) = 0. 

The time-independent creation and annihilation operators obey the canonical anticommutation rules 

Ka'^a'} = 2 60 ( 2?r ) 3 ^ Sx A ' ' 

{djf,A»4>,A'} =2e ° S $-$) ■ 

Following the method of refs. [ll|, G3 , it proves convenient to write 

-. _i r n -> i 

Ux(k,rj) = (e + ma dc ) 2 d v - 7 • k + m a(rj) f k (rf)Ux 

V\(k,rj) = (e + ma dc ) 2 87 d n + 7 • A; + m a(rj) Va , 

(eo + to adc) 1 being a normalization factor and Ux, Va being constant spinors [III ll2T| obeying 

fUx=Ux , 7°V A = -V A , A = ±l. 

More explicitly, 



U\(k,rj) = (e + m a dc ) 1 
Vx(— fc, r?) = (e +m a d cY 



[i f'M + m a( V ) f k {rf)} 

x k f k {n) 



Ux , 



\kg k (j]) 

H.9fc( ? /) +" la (l) 9k(v)] 



(3.4) 



(3.5) 



(3.6) 
(3.7) 



(3.8) 

(3.9) 
(3.10) 

(3.11) 



(3.12) 
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The mode functions fk(v)i 9k{?l) obey then the following equations of motion 

drj 2 

- £_ 

dvj 2 



„ + fc 2 + m 2 a 2 (77) — i m a!{rj) 



+ k 2 + m 2 a 2 (rj) + im &'(?7) 



A fa) 

Sfefa) 



(3.13) 
(3.14) 



We choose the initial state for the fermion field as the vacuum state which is a thermal equilibrium state at temperature 
T for the fermion. This Fock state |0 > is annihilated by the operators 6^ , and , . 
Therefore, we have as initial conditions for the mode functions ujl LL2I 



A(0) = 1 , f' k (Q) = -ie , 
fl fc (0) = l , g' k (0) = +ie 



(3.15) 



These initial conditions describe the Bunch-Davies vacuum when they are applied at asymptotically earlier times in 
the past (ij — > — oo) @, See the discussion in sec. I VIII below. 
Eqs.([3 ^ -(|3T5 ]) imply that 

9k(rf) = fk(v) ■ 

That is, we have only one independent and complex mode function. 
The scalar products of the spinors U\(k, 77), V\(k, rj) take the values 



Ul(k,J]) Uy(k,7 ] )=2e S x v 

v£(k,ri)V x ,(k,ri)=2e S x y 

As a consequence, the mode functions obey the relation [HI \l2^ 

\f k ( V )\ 2 - i m a( V ) [f k ( V ) f'Xrj) - f k (v) fUv)} + [k 2 + m 2 a 2 (77)] \f k (77)! 2 = 2 e (e + m a dc ) 

which provides a conserved quantity. 

The energy momentum tensor for a spin 1/2 field is given by Q 



T, 



^ 2 

and its expectation value has the fluid form 

< T F ° > = < H F > (77) , < T F l > = -Sj <P F > (77) 

since we consider homogeneous and isotropic quantum states and density matrices. 
More explicitly, the energy density in conformal time takes the form 

< H F > (v) = < *(^,?/) f H F ®(x,rj) > , 
where the fermion hamiltonian Hp is defined by 

a(rj) 70 Hp = —tf ■ V + m a(rf) = 7 • p + m 0(77) . 
An analogous expression can be written for the pressure, 

1 



<P F > (77) 



< * 7 • p * > (77) 



3 a(r)) 

Here too, it is convenient to consider the conformal energy and pressure, 

e F {ri) ee a 4 (77) < Hp > (77) , 7^(77) = a 4 (77) < Pp > (77) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



(3.21) 
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We find the trace of the energy-momentum tensor from eqs. (|3.19p . (|3.20|) and (|3.2ip . 

e F (rj) - 3 pf(v) = m a{rj) Sf(??) , or (3.22) 

< H F > (?y) - 3 <P F > (77) = to < > (77) . 
This is the expression of the virial theorem in the present context and 

E F (»j) = a 3 ( V ) < > fa) . (3.23) 
The above expressions for the energy density and pressure obey the usual continuity equation in cosmic time 

^<H F > +3 H(rj) (< H F > + < Pf >) = , (3.24) 

In conformal time by using eqs. (|3.22p - (13.23[) the continuity equation (|3.24p becomes, 

de F da(rj) 

— — = m — — T, F (r]) . (3.25) 
dr\ arj 

We thus see from eqs. (|3.22p and (|3.25p that there is only one independent quantity among s f (t]), V F {rf) and Y, F (rj). 

IV. THE COSMOLOGICAL QUANTUM VACUUM. 

There are two widely separate scales in the field evolution in the cosmological space-time. The fast scale is the 
microscopic quantum evolution scale, typically ~ 1/M ~ 1/m. The slow scale is the Hubble scale 1/H of the universe 
expansion. 

When M ~ m 3> H we can consider that the scale factor is practically constant. Therefore, in conformal time the 
quantum field evolution is like the evolution in Minkowski space-time with a mass M a(rj) or m a(rj) for bosons or 
fermions, respectively [see ecis. (|2.5p and (|3.4[) ] . 

The scalar and fermion densities follow as equal point limits of the scalar and fermion t wo p oint functions. That 
is, we consider the scale factor a as a constant and obtain for the scalar two point function [14j |. 

G s (x-x' )V - V ',Ma) = <T<p(2,r,)<p& t tf)> = ± ! f ^ «-*''-»'■ 

1 M 



a 2 J (2tt) 4 k 2 -a 2 M 2 + 10 



— K x {Maz) , zee v /(f_f')2_( ?7 _ ?7 ')2 ) 



(2tt) 2 za 

where Kx (x) is a modified Bessel function [32j ■ 

Eq. (|4.ip is the zeroth order adiabatic approximation. It differs from the exact two point function by quantities of 
the order 0(a'(T])), 0(a"(rj)) etc. 

We find from eq. (|4.1|) in the coincidence point limit [32| : 



log (M a z) + C - In 2 - - 



[1 + 0(M 2 z 2 )] . (4.2) 



where C = 0.57721566 ... is the Euler-Mascheroni constant. Eqs. (|4.1| )- ()4.2|) display the two point functions for the 
zero temperature vacuum. The effect of a non-zero temperature on the two point function is negligible for a> 1 as 
we show below [eg. (|4.10|> ] . 



The fermion two point function takes the form 14 1 



and hence, 

G F (x — £ , rj — 77 , ma) = < T ^(x, if) ty(x , r\ )> = —4m Gs(% — x , r\ — 77 , ma) , Dirac Fermions . (4-4) 
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The minus sign in front arose from the anticommutation of the fermion fields going from eq. (|4.3[) to eq. (|4.4p . Here we 
used ea. (|4.1[) and 



Tr U = 



Tr 1 = 4 



(4.5) 



That is, the factor 4 = 2 x 2 in eqs. (|4.4p and (|4. 5|) comes from the fermion and antifermion contributions times the 
number of helicities of a Dirac fermion. Hence, this factor 4 becomes a factor 2 for Majorana fermions: 

Gp(x — x',r) — rj',m a) = < ^(x,rj) ^(x',rj') > = —2 m Gs(x — x! ,rj — rj ,ma) , Majorana Fermions . (4.6) 

We find in the coincidence point limit [32| : 



n ,~ -n i \ 2 K to 

Gf{x — x ,?7 — rj ,ma) = — - 



1 



1 



log (to a z) + C — In 2 — 



1 



(2tt)2 

Here, H = 1 for Majorana fermions and H = 2 for Dirac fermions. 
In order to define the vacuum densities as the coincidence limits, 

< if 2 > {rj) = < V 2 (f,7?) > , < > (77) = < ^>(X, rj) ^{x,i]) > 



[1 + 0(m 2 z 2 )] . (4.7) 



we have to subtract the singularities at z — in eqs. (!4.2|) and (|4.7| . Subtracting the singularities leaves a finite z 
independent piece. Requiring that the vacuum densities vanish in Minkowski space-time (a — I) we obtain, 



<<p 2 > (77) 



M 2 H to 3 

^^[loga + bsMa)} , <**>(„)=-_ 

The functions fs(o-) and /f(o) are finite and vanish for Minkowski space-time, 

/s(i) = o , Mi) = o. 

We compute in [3§| the terms bs fs{ a ) an d bp fp(a) with the result: 



[log a + b F f F (a) 



(4.8) 



/s(oo) = 1 + 



j. , > a(r))>a dcB ,a dcf 

, /f(oo) = 1 + 



When one performs an infinite subtraction at z = 0, an additional finite subtraction can always be done. We recognize 
that the additional terms containing 65 and bp can be absorbed in a finite multiplicative renormalization of the scale 
factor. That is, introducing 65 and bp amounts to a scale transformation. We compute in [3^| the coefficients 65 and 
bp in terms of the subtraction scale in momentum space (x M) for scalars and {x m) for fermions, with the result 



b s {x) = b F (x) 



1 



log x - log a dc 



where a^c stands for the scale factor at decoupling time (initial time). In summary, we have for the late time regime, 



M 



. 2 - / \ a(r;)>a dcs ,a dcf M 2 . 

<ip>M = w*? [ga{n)+s] ~n^? 

<W> ( V ) ■<"»*"— [loga W + 6f] N 



log 



a(ry) 



1 



(2^ 



(2 7T) 



log 



o(r?) 1 
x a dc f 2 



(4.9) 



where adcs and aa c f stand for the scale factor at the decoupling times (initial times) for the scalar and the fermion, 
respectively. 

The two-point functions eqs. (|4.ip and (|4.4p correspond to the zero temperature case. The singular pieces for z — > 
are temperature independent. We can disregard the temperature dependent contributions to the two point functions 
since for large a they decrease as 



'M 



T 
2 7r a 



a > 1 



(4.10) 



12 



The scalar and fermion densities < ip 2 > (77) and < 'J'S > (77) can be also computed as momentum integrals over 
the mode functions 4>k(v) an d fk(v)- This computation is done for late times in a forthcoming publication (39j and 
confirms the results presented here. In addition, we obtain the subdominant terms in l/a 2 (r]), a(rf) / 'a 2 '{rf] , . . . etc. 

The equal points behaviour of the two point function eqs. (|4.2[) and (14. 7|) is generic for any curved space-time when 
expressed as a function of the geodesic (squared) distance a = z 2 a 2 between the two points. That is, the short 
distance behaviour is uniquely and universally determined by the local space-time geometry. It must be noticed that 
the divergences and finite pieces at a = are of the same type as in Minkowski space-time. This is the so called 
Hadamard expansion for a — > and is equivalent to the adiabatic expansion. The coefficients of the divergent and 
finite parts are called Hadamard coefficients and they are known for generic space-times. This will be discussed in 
detail in connection with the cosmological dark energy evaluation in ref. [39j . 



A. Vacuum Energy Density and Pressure for Late Times 

The total energy density e{ri) and pressure V{rf) 

<H>{rj) = <H s >{ri) + <n F >{rj) , < P > (77) = < P s > (77) + < P F > (v) (4.11) 



can be computed in the late time regime using the virial theorem eqs. (|2.15|) and (|3.22|) . the continuity equation eqs. 
(|2.17p and (13.25|) and the late time behaviour of the densities, eq. (|4.9p . 



We obtain after calculation for the energy density and pressure, 



< n > (v) 



<p>(v) 



M 4 



a(>j)3>a dcs ,a dc f 



2(4tt) 2 
M 4 



log 0(77) +b s 



1 



(4tt) 2 



2(4tt) 2 



loga(rj) +b s + 



12 



(4tt) 2 



log a(7?) + bp 



loga(ry) +b F + 



12 



(4.12) 



The decoupling (initial) times for the evolution of scalars and fermions can be different to each other. We have 
absorbed in bs and b F the corresponding initial values of the scale factor for scalars and fermions, respectively. 
The positivity of the energy density impose the condition, 



M 4 >2Km 4 



Notice that, 



< p > (77) + < n > (77) 



a(i))S>a dca ,a dcf 



6(4tt) 2 



[M 4 — 2 H m 4 ] 



is time independent and independent of the finite subtraction coefficients bs and b F as well. 
From eq. (|4.12p we obtain for the equation of state, 

, s < P > (?7) a( t? )»a dcs ,a dct 1 1 



log a (77) - 



1 



v 2Km 4 l 
OS j-pr- Of 



1 



2 Km 4 
M 4 



That is, we find w(ri) < — 1 with 711(77) asymptotically reaching the value —1 from below. 
It is convenient to express the scale factor in terms of the redshift as 



0(77) e 6 



1 + z 



a(rj) e bF = 



l + z F 
1 + z 



(4.13) 



(4.14) 



where zs (z F ) is the redshift when the evolution of the scalar (fermion) become the one of a free field in the cosmological 
space-time. In terms of z$ and z F eas. (|4.12l) read, 



< H > (z) 



2(4tt) 



M 4 log(l + z s ) - 2 H m 4 log(l + z F ) — (M 4 — 2 H m 4 ) 



log(l + z) 



1 



(4.15) 
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< P > (z) = 



M 4 log(l + zs) — 2 H to 4 log(l + z F ) - (M 4 - 2 H m 4 ) 



2(4tt) 2 

The equation of state (|4.13[) as a function of z takes the form 



log(l + z) - 



12 



w(z) + 1 



1 - 



2_N_m_ 

M 4 



3 log(l + * s ) - ^ log(l + zf) - (1 - [log(l + z) + i] 

The equation of state and the energy density become today: 



w (z = 0) + 1 = 
< W > (z = 0) 



1 



2 K m 
M 4 



3 log(l + z<?) 4 
1 



1 



2(4?r) 



Af 4 



log(l + z s ) 



2 H m 4 



log(l + z F ) 



The energy density at late times r\ after decoupling and the energy density today are related by 

rti ( \ a(r()>a d c 3 ,adcf oy ^ / Af 4 - 2 K m 4 0(77) 

<H> (77) = < « > (2 = 0) H ^ — log ■ 



2(4tt) 2 



«o 



.(4.16) 



(4.17) 



(4.18) 



(4.19) 



where we used eas. (|4.12p and (|4. 1 8[) and ao stands for the scale factor today. We identify the vacuum energy density 
today < TC > (z — 0) with the observed dark energy p\. We can then write, 



< H > (??) = pa 



1 + (3* loj 



a(r?) 



a 



(4.20) 



where 



M 4 



log(l +z s )-\- ^ [log(l + z F ) - i] 



(4.21) 



That is, the vacuum energy density at late times after decoupling grows as the logarithm of the scale factor. Moreover, 
the equation of state approaches — 1 from below as, 



w{rj) + 1 



a(jj)>a d cs,adcf 



M 4 - 2 H TO 4 



6(4^) 2 p A 



1 + A» lot 



a ('7) 



The previous equations in this subsection generalize when there are several scalar and fermion fields just summing over 
their respective contributions. Let us consider the case of several scalars and fermions. This case is relevant to study 
whether the three neutrino mass eigenstates can contribute to the dark energy. Eqs. (|4.18[) become for zs, zp 3> 1, 



Y.4 M t - 2 ^ E 



4 



Pa 



2(4tt) 2 



6 (4ir) 2 p\ 

logz s -£)5> 4 -2K (log z F -£)]T, 



(4.22) 
(4.23) 



where j and i label the species of scalars and fermions, respectively. It is convenient to eliminate the sum of scalar 
masses ^ ■ M 4 between eas. (|4.22l) and (|4.23[) with the result, 



w(z = 0) + 1 = 



log z s 



ZF 



3(4^) 2 pa 



(4.24) 



We see here that the scalar contributes to the equation of state today by the negative term — 1/[3 (log zs ~ j)] while 
the fermions give for zs > zf a positive contribution proportional to the sum of the fourth power of their masses. 
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B. The Quantum Nature of the Cosmological Vacuum. 

Local observables as < (p 2 >, < 'i'vp >, the energy density and the pressure involve the product of the field operators 
at equal points. This is identical to one-loop tadpole Feynman diagrams. Logarithmic dependence on the scale of 
the momenta is typical in one-loop renormalized Feynman diagrams [l4|]. Here, we analogously find a logarithm of 
the scale factor in eqs. (|4.9p and (|4.12[) through the same mechanisms at work in renormalized quantum field theory. 
Hence, the dark energy follows here as a truly quantum field vacuum effect. We stress quantum field effect and not 
just quantum effect because the infinite number of filled momentum modes in the vacuum as well as the subtraction 
of UV divergences play a crucial role in the vacuum late time behaviour. Here the quantum fields are not coupled 
neither self-coupled but they interact with the expanding space-time geometry. 

Notice that these results eqs. (|4.9|) . (|4.12p and (|4. 13|) are valid for any expanding universe. They do not depend on 
the specific time dependence of the scale factor 0,(77) provided it grows with 77. 



The quantum nature of the these vacuum cosmological effects in the physical observables are manifest from eqs. (|4.9p 
and (|4~T2| . 

M 2 c 2 M c 

<<p 2 >(ri) ~ M 2 log 0(77) = — — - log 0(77) = — — log 0(77) , 

h A c 

-j^J log 0.(77) = -jf- loga(?7) . (4.25) 

These quantities are of quantum nature since they depend on h. There is no 'classical contribution' to the vacuum 
energy. Eq. (|4.25[) just means that the scale of the dark energy density is of one scalar rest mass per a volume equal 
to the cube of the Compton wavelength Ac for the scalar particle. Notice that Ac = -A— — 0.05 mm is almost a 
macroscopic length while the mass of the scalar particle M ~ 4 meV = 7.1 10~ 36 g is extremely small (see below for 
the value of M). 



V. DARK ENERGY FROM THE COSMOLOGICAL QUANTUM VACUUM. 

The current value for the dark energy density is [2q . 

Pa = n APc = 3.28 x KT 11 (eV) 4 = (2.39 meV) 4 , (5.1) 

corresponding to h = 0.73 and J7a = 0.76. Here 1 meV = 10~ 3 eV. 

Bosons give a positive contribution to the dark energy through the cosmological quantum vacuum while fermions 
give a negative contribution. Therefore, the boson contribution must dominate. 

As discussed in sec. IV Al the lightest neutrino certainly contributes to the cosmological quantum vacuum unless it 
dissipates. Definitely, a boson contribution is needed. The photon and graviton contributions are irrelevant since 
their masses are most probably zero and at most [25[ 

TO 7 < 6 X 10" 17 eV , m gravi ton < 7 x 10~ 32 eV . 

Massless particles contribute to the energy-momentum tensor through the trace anomaly @, Q • This contribution is 
of the order of Hq where Hq is the Hubble parameter today. Its value 

H = 1.558 x 1(T 33 eV , (5.2) 

is exceedingly small to explain the observed value of the dark energy. 

A scalar particle can produce the dark energy eq. (|5.ip through its quantum cosmological vacuum provided: 

• Its mass is of the order of 1 meV and it is very weakly coupled. 

• Its lifetime is of the order of the age of the universe. 

Spontaneous symmetry breaking of continuous symmetries produces massless scalars as Goldstone bosons. If, in 
addition, this continuous symmetry is slightly violated the Goldstone boson acquires a small mass. This is the natural 
mechanism that generates light scalars and several particles have been proposed on these grounds in the past. The 
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axion is certainly the one that caught more attention in the literature. Other proposed particles are the familons and 
the majorons [2ll.[22l]. 

The (invisible) axion [l7| (if it exists) is hence a candidate to be the source of dark energy. 

Axions were proposed to solve the strong CP problem in QCD [16j |. Axions acquire a mass after the breaking of 
the Peccei-Quinn (PQ) symmetry when the temperature of the universe was at the PQ symmetry breaking scale ~ f a 
fll ]. All axion couplings are inversely proportional to f a and the axion mass is given by [l8j |. 

„ /10 9 GeV\ „ ,„ nS 
M a ~ 6 x f J mcV . (5.3) 

The following range ('axion window') is currently acceptable for the axion mass (l9j : 

1CT 3 mcV < M a < 10 meV . (5.4) 

Therefore, this pseudoscalar particle has extremely weak coupling to gluons and quarks and hence it contributes to 
the cosmological quantum vacuum. For example, the axion-photon-photon coupling is given by [H, [HI] 



9a 7 7 



I0- 10 / 



GeV V 1 mcV / 



As a consequence, the axion lifetime to decay into photons is much longer than the age of the universe. Dissipation 
of the energy in the cosmological quantum axion vacuum takes longer than the age of the universe too. 

An axion with mass ~ 1 meV and hence f a ~ 10 9 GeV decoupled from the plasma at a scale of energies ~ 2 x 10 5 GeV 
[20l |. that is z$ ~ 2.2 x 10 18 . The temperatures of the axions and neutrinos today are lower than that of photons 
today, 

TUoday = ^yj) Tcmb today = 0.1676 meV , T a today = 0.078 meV . (5.5) 

Since the axion lifetime is of the order or larger than the age of the universe, no specific properties of the axion play a 
role in the dark energy except its mass and decoupling redshift. However, the dark energy depends on the decoupling 
redshift rather weakly since it is through the logarithm [see eq. (|4.15p ]. 

Neutrinos in the universe are believed to be effectively free particles when the temperature of the universe is below 
~ 1 MeV @. That is, neutrinos decouple at a redshift zp ~ 0.6 x 10 10 . Before, such time electrons and neutrinos 
interact, keeping them in thermal equilibrium. 

Therefore, we can treat the axion with mass ~ 1 meV and the lightest neutrino as free particles in the universe for 
redshifts z < zg ~ 2.2 x 10 18 and z < z F ~ 0.6 x 10 10 , respectively. 



A. Neutrino Mass Eigenstates 

As is known, the two heavier neutrino mass eigenstates v-i and 1/3 with masses mi and 7713, respectively, annihilate 
with their respective anti-neutrinos yielding the lightest neutrino eigenstate v\ and its antiparticle through weak 
interactions. However, this process is too slow for nonrelativistic neutrinos even compared with the age of the 
universe. Their decay rates can be estimated to be 

r, ~ Gl ml \~ — , r 3 ~G 2 P r4- ' 



1.5 x 10 33 yr ' f J 5 x 10 29 yr ' 

where Gf = 1.166 1 - 23 (eV) -2 stands for the Fermi coupling. 

Neutrinos with masses ui2 ~ 0.01 eV or 7713 ~ 0.05 eV will produce through their cosmological quantum vacuum 
today a large negative contribution to the dark energy. 

Therefore, the heavier neutrinos {v 2 and ^ 3 ) must annihilate with their respective anti-neutrinos into the lightest 
neutrino v\ through a mechanism such that 

r 3 > T2 > age of the universe -1 . (5-6) 

Our results for the dark energy are independent of the details of the decay mechanism. All what counts is that the 
decay rates of the heavier neutrinos obey eq. (|5.6p . 
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As a minimal assumption, let us consider the following effective couplings between the neutrinos, 

— ^ 2 1-2 *i *i , ^3 *3 *i *i , (5.7) 

where M' is a mass scale much larger than the neutrino masses. We thus find, 



2 ~ M / 4 ' 3 ~ M /4 



Imposing eq. (|5.6|) yields, 

M' < 1 MeV , for m 2 = 0.01 eV and M' < 10 MeV , for m 3 = 0.05 eV . (5.8) 

The first estimated bound (1 MeV) applies for a direct hierarchy of neutrino masses (777,3 ~ 0.05 eV > 7772 ~ 0.01 eV > 
mi ) while the second estimate (10 MeV) is for an inverse hierarchy of neutrino masses (7773 ~ 7772 ~ 0.05 eV > mi). 

Effective couplings of the type eq. (|5.7[) can be obtained from different renormalizable models. 
Notice that the two heavier neutrinos decays contribute to the background of lighter neutrino particles and not to 
the neutrino quantum vacuum. 

Lagrangians leading to effective couplings analogous to eq. (|5.7|) have been considered in the context of models to 
generate neutrino masses and to provide light dark matter candidates [36] . Moreover, mass ranges compatible with 
eq. (|5.8p have been obtained from various and independent considerations [H, [36J. In case the effective couplings 
eq. (|5.7[) arise from Yukawa couplings of the neutrinos with a scalar particle of mass M', this scalar particle cannot be 
a dark matter candidate since it decays into neutrino- antineutrino pairs. 

The lightest neutrino with mass mi can be self-coupled through the interaction 



1 



M" 2 

Its decay rate, 

Ti 



M" 4 

is of the order or larger than the age of the universe when 



5 



M" <( — Y 50keV. (5.9) 
V meV / 

Hence, if eq. (|5.9p is fulfilled, the energy in the neutrino vacuum dissipates into the lightest neutrinos v\, thus con- 
tributing to the neutrino background. 



B. Light particle masses and the dark energy density today 

Let us consider the case where only one light scalar field and one light fermion field contribute to the quantum 
vacuum energy. That is, a light scalar and the lightest neutrino. We obtain from eq. (|4.18|) for the mass of the scalar, 



M 



2* V^pI 

(\ogz s -\) 



H 777 4 



(4 tt) 2 p A 



\ogz F 



(5.10) 



where we identified the vacuum energy density today < TL > (0) with the observed dark energy p\. 

We now obtain using the observed value of the dark energy eq. (|5.1[) and the decoupling redshift for the neutrino 
z F - 0.6 x 10 10 , 



M 



10.1 meV 

{\ogzs-\f 



1 + N 



/ to \ ' 
V3.90meV/ 



(5.11) 
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If the lightest neutrino has a very small mass [29| m <C 1 meV or if it decays in the time scale of the age of the 
universe [see eq. (|5.9p ] so the neutrino vacuum dissipates, there is no neutrino contribution to the dark energy. In 
these cases eq. (|5.11|) gives for the mass of the scalar: 

10.1 meV . , „, 

M = , no vacuum neutrino energy . (5-12) 



Assuming the scalar field to be the axion we can use the value zs ~ 2.2 x 10 18 for the axion decoupling redshift 
and eq. (|5.1ip becomes, 



M(m) = 3.96 meV 



1 + 



3.90 meV 



(5.13) 



The values of the neutrino masses are not yet known, only their differences arc experimentally constrained [22j. Both 
in the direct and inverse mass hierarchies the mass m of the lightest neutrino can be in the meV range (or even to be 
zero) . 

According to ref. [26}, we have 

1 

m = — m,2 

O 

where rri2 is the mass of the middle neutrino. Combining this with the known neutrino mass differences [22j yields 

M 

m = 3.2 ±0.1 meV (5.14) 
This value for the neutrino mass perfectly agrees in order of magnitude with the see-saw prediction [22|. 

M 2 

- Formi ~ 6 x 10~ 3 eV , 
Mqut 

for the typical values Afpcrmi = 250 GeV and Mqut = 10 16 GeV, of the Fermi and Grand Unified energy scales, 
respectively. 

Eqs. (|5.13[) - (|5.14p gives for the axion mass: 

M(m = 3.2 meV, N = 1) = 4.35 meV , M(m = 3.2 meV, K = 2) = 4.66 meV , (5.15) 

for Majorana and Dirac neutrinos, respectively. 

If the lightest neutrino has a very small mass[Ii| mi <C 1 meV or if it decays in the time scale of the age of the 
universe [see eq. (|5.9|) ] so the neutrino vacuum dissipates, there is no neutrino contribution to the dark energy. In 
those cases the axion mass is given by 

M = 3.96 meV , no vacuum neutrino energy . (5.16) 

All the axion mass values eq. (|5.13p and eqs. (|5.15p - (|5.16[) describe the dark energy observed today eq. (|5.ip . The 
numerical values for the axion mass in eqs. (|5.15p - (|5.16p are within the astrophysical bound of eq. (|5.4p . 

We compute the equation of state today from eq. (|4.18p and display it in table 1 in three relevant cases: (i) no 
neutrino contribution to the dark energy, (ii) a Majorana neutrino contribution, (hi) a Dirac neutrino contribution. 
In all three cases the observed value eq. (|5.ip of the dark energy is imposed. For the last two cases we choose for the 
neutrino mass m = 3.2 meV and the scalar mass given by eq. (|5.15p . 

We see that w(0) is slightly below —1 by an amount ranging from —1.5 x 10~ 3 to —8 x 10~ 3 . 

It can be noticed that the mass of the lightest neutrino [eq. (|5.14p ] turns to be much higher than today's neutrino 
temperature: 

m Dirac -, n r m Majorana no „ /_ 1 _> 

= = 19.6 , — H = 23.3 , (5.17) 

^ ftoday -*■ ftoday 

where we used eg. (I5.5[) . That is, the neutrinos forming the neutrino background are today non-relativistic particles. 
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Neutrino Type 


Scalar Mass 


Equation of state today 


no vacuum neutrino energy 


M = 3.96 meV 


to(0) + 1 = -.00794 


Majorana Neutrino 
m = 3.2 meV 


M = 4.35 meV 


w(0) + 1 = -.00473 


Dirac Neutrino 
to = 3.2 meV 


M = 4.66 meV 


w{0) + 1 = -.00156 



TABLE 1. The equation of state today w(0) + 1 computed from eq. ()4.18p in three relevant cases which all describe 
the dark energy observed today [eq. (|5.ip ]: (i) no neutrino contribution to the dark energy, (ii) a Majorana neutrino 
with mass m = 3.2 meV, (hi) a Dirac neutrino with mass to = 3.2 meV. See the discussion in sec. IV Bl 



Let us now analyze the possibility that all three neutrino eigenstates contribute to the dark energy. This contribution 
crucially depends on the values of their masses to the power four through the dimensionless factor 



1 Ei 



in 



4 



3(4^) 2 p A 



as we see from cq. p~2"2>P~2"4"|) . 

For the normal hierarchy we have [2z 



mi = 3.2 meV , TO2 = 9.5 meV , to 3 = 47 meV , 
and for the inverted hierarchy [22j 

nil = 3.2 meV , TO2 = 47 meV , 7713 = 48 meV . 
Thus, using eqs. (|5.ip the factor T takes the values 

J : normal = 315 , ^inverted = 656 . 

Inserting these numbers in the equation of state today eq. (|4.24p yields values for w(0) in strong disagreement with 
the data unless we fine tune z$ — zp. Since there is no reason to have such equality, we conclude that the vacuum of 
the two heavier neutrinos must not contribute to the dark energy. Their quantum vacuum must dissipate as discussed 
in sec. IV Al 



VI. THE FUTURE EVOLUTION OF THE UNIVERSE 



The future evolution of the universe follows by inserting the total energy density in the Einstein- Friedmann equation 



Ht , 



where we use cosmic time t, G is the gravitational constant and the total energy density H.t is the sum of the 
contributions from the dark energy, the matter and the radiation. 

We obtain using the dark energy expression eq. (|4.20p the self-consistent Einstein- Friedmann evolution equation, 



H\t)=Hl 



n A 1 + /?n ioi 



,o(t) 
a 



a, 



3(t) 



ad 



(6.1) 



where a = a (today), is defined by ea. (|4.2ip . p A = p CIlt Q A is given by eq.flTS]), being p crit = |^%, H 
h/ [9.77813 Gyr] stands for the Hubble parameter today, £l A — 0.76 = 1 — fimatter — ^rad |25j. 
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We get using the explicit values for M and m eqs. (15.14l) - (|5.16l) : 
Po = 0.0238 no vacuum neutrino energy , j3\ = 0.0347 Majorana neutrino , 02 — 0.0459 Dirac neutrino 
For a(t) > ap, the matter and radiation contributions can be neglected in eq. (|6.1|) and we have, 



d log a(t) 



dt 

This equation can be immediately integrated with the solution 

Ho t>i 



i + A<log^ 



a(t) a e Cl H ° t+C2 {H ° 4) . (6.2) 



where 



ci = = 0.87 , c 2 = i fl A fa = 0.19 fa , 

0.00452 < c 2 < 0.00872 . (6.3) 

The left and right ends of the interval in ci eq. (|6.3p correspond to the cases in which there is no neutrino contribution 
and to the lightest neutrino being a Dirac fermion with mass m = 3.2 meV, respectively. 

We find that the universe is presently reaching an asymptotic phase where it expands as indicated by eq. (|6.2|) . 

Eq. (|6.2[) shows that the expansion of universe is in the future faster than in the de Sitter universe. 

Notice that the time scale of the accelerated expansion is huge ~ 1/Hq — 13.4 Gyr. The quadratic term dominates 
over the linear term in the exponent of ea. (|6.2|) by a time t ~ 100/i?o to 200/i?o- 

In this accelerated universe, eq. (|6.ip shows that the Hubble radius 1/H decreases with time as 

1 1 



H H v/Togo(i) 
VII. DISCUSSION 

The non-trivial energy and pressure that we have is an effect resulting of the expansion of the space-time as it arises 
from the log 0(77) factor in eqs. (|4.12p . No dark energy appears in Minkowski space-time. Namely, the formation and 
growth of the vacuum density, the vacuum energy density and pressure is an effect due to the presence of quantum 
fields in an expanding cosmological space-time. 

Notice that the energy scale of the cosmological vacuum is given by the mass of the particle when this mass is larger 
than the Hubble constant [see eq. (|5.2[) ]. For massless particles, the energy scale of the cosmological vacuum is given 
by the Hubble constant. 



The axion evolution for z > 10 18 as well as the neutrino evolution for z > 10 10 are beyond the scope of this article. 
Namely, the regime where the interaction of axions and neutrinos with the plasma particles cannot be neglected. 
We choose as initial state for both the axions and the neutrinos the vacuum thermal equilibrium state. It must be 
remarked that the vacuum energy at late times is independent of the initial temperature as shown by eq. (|4. 10[) . 

Before decoupling, particle interaction is non-negligeable and dissipation is important depleting the vacuum energy 
,Hlj|. Hence, the vacuum energy can only become significant after decoupling. Therefore, it is a good approximation 
to just study the free quantum field evolution in the cosmological space-time after decoupling. 

The initial conditions eqs. (|2.9p and (|3.15[) are imposed at the origin of the conformal time. We shall see now that 
they are equivalent to the Bunch-Davies vacuum conditions. Since the initial time corresponds to a large value of 
redshift, it corresponds to asymptotic times in the past in a very good approximation. More precisely, the conformal 
time is related with the redshift z by 



?y = —j=^ , matter dominated era , 
V 1 + z 
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2t ^TTl- t radiation dominated era (7.!) 

l + Z yTT^ 

where to = 13.7 Gyr is the age of the universe, 1 + z cq = 3048 is the transition from the radiation dominated to the 
matter dominated era. 770 = 3 to corresponds to the present time. For z z cq we see that, 

= 0-018 t . 

V 1 + z cq 

Hence, the conformal time at decoupling differs from the conformal time today 770 = 3 to by an amount ~ 3 to- Hence, 
the initial time can be considered as an asymptotic time deep in the past. More precisely the change in the phases of 
the mode functions is characterized by M to ~ 3 x 10 30 for a typical mass M ~ 4 meV. Hence, the initial conditions 
for the mode functions eqs. (|2.9[) and (|3.15[) are virtually identical to Bunch-Davies initial conditions. 

The vacuum density and energy density eqs. (|4.9p and (|4.12|) are determined by the short distance behaviour of 
the two point function in coordinate space. In momentum space, it is the high energy behaviour that dominates the 
vacuum density and energy density for late times. The physical quantities can be written as integrals of mode functions 
as in eqs. (l2.7p and (|3.5p . One can see that the relevant comoving momenta k's contributing at a physical energy scale 
q take the value k = q a{rf). At late times, (e. g. today) a(rj) ~ ^decoupling and so only large k ~ ^decoupling M are 
relevant [3!| . This fact suppresses the effect of the initial conditions on physical quantities measured today. Analogous 
effects take place for the initial conditions of inflationary fluctuations (with the exception of the quadrupole) [24[ . 
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FIG. 1: The equation of state w(z)vs.z for the three cases explicitly calculated in this paper: (i) [full line] no neutrino 
contribution to the dark energy and the scalar mass M = 3.96 meV, (ii) [broken line] a Majorana neutrino with mass m — 3.2 
meV and the scalar mass M = 4.35 meV, (iii) [dotted line] a Dirac neutrino with mass m = 3.2 meV and the scalar mass 
M = 4.66 meV. [See the discussion in sec. IV B] , In all three cases w < — 1 by less than 1%. 



In fig. [T]we plot the equation of state w(z) as a function of z for the three cases explicitly calculated in this paper: 
(i) no neutrino contribution to the dark energy and the scalar mass M — 3.96 meV, (ii) a Majorana neutrino with 
mass m — 3.2 meV and the scalar mass M — 4.35 meV, (iii) a Dirac neutrino with mass m = 3.2 meV and the scalar 
mass M = 4.66 meV. [See the discussion in sec. IV B| . 

We see that the equation of state in all three cases differs from the cosmological constant case w — — 1 by less than 
1%. 
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Several recent papers claim that values of the equation state w a little below —1 are favoured by the observations 
[23l ] . This is precisely like the equation of state eq. (|4.17p found here. However, the observed values w < — 1 are within 
the error bars and so the present data do not give a definite conclusion on whether w < — 1 (33l |. 



The value of the lightest neutrino mass eq. (|5.14|) is well below the neutrino mass splittings y/ Am^ un and \J Am^ tm 
and consistent with both direct and inverse mass hierarchies [25[ . A quasi-degenerate mass spectrum will give a large 
negative contribution to the dark energy and will require a scalar particle with a mass M > 100 meV to reproduce 
the observed dark energy data ea . (| 5 . 1 1) . Such a particle can very well exist but it cannot be the axion [see ea. Q5.4p ]. 
Indeed, the scalar particle can have the mass value given by ea. Q5.16[) in case all three neutrinos decay in a time scale 
of the age of the universe in order to dissipate their cosmological vacuum energy as discussed in sec. IV Al 

It is claimed in ref.[3l| that successful baryogenesis is possible for a range of neutrino masses from 10~ 3 eV to 0.1 
eV in agreement with neutrino mass differences from oscillations and the value eq. (|5.14p for the mass of the lightest 
neutrino. 

A weakly coupled scalar particle with a mass in the range 1 — 1.5 meV is proposed to explain the PVLAS experiment 
[3^ . This mass is, however, 3 — 4 times too low to produce the observed dark energy. 

VIII. CONCLUSIONS 

We find that the presence of a cosmological quantum vacuum energy with an equation of state just below — 1 
is the unavoidable consequence of the existence of light particles with very weak couplings. Bosons yield positive 
contributions and fermions yield negative contributions to the vacuum energy. 

It must be noticed the present lack of knowledge about the low energy (energy ~ 1 meV) particle physics region. 
Actually, most of the constraints on this sector follow from astrophysics and cosmology [18|]-[2(| including the new 
constraints that we obtain here on the axion mass. 

No exotic physics needs to be invoked to explain the dark energy. Since the observed energy scale of the dark energy 
is very low, we find natural to explain it only through low energy physics. The effects from energy scales higher than 
1 eV or even 1 MeV arrive strongly suppressed to the dark energy scale of 1 meV. 

In summary, dark energy can be explained by a very light and very weakly coupled scalar particle which decouples 
by redshift z$ 3> 1. If the scalar particle is the axion, then z$ ~ 2.2 x 10 18 . 
We have four main cases: 

• (i) No neutrino contribution. This happens when the lightest neutrino has a mass m< 1 meV and when the 
vacuum neutrino contribution dissipates in the time scale of the age of the universe [see ea. Q5.9p ]. The scalar 
mass must be 

10.1 meV . ,„„, 

M = , no vacuum neutrino energy . (8-1) 

(l0g25-!) 3 

If the scalar is the axion, then M = 3.96 meV in this case. 

• (ii) The lightest neutrino is Majorana and has a mass m ~ 3.2 meV. Then, the scalar mass must be 

11.1 meV , ™ . 

M = 5- , the Majorana neutrino contributes. 

(logzs-lf 

If the scalar is the axion, then M = 4.35 meV in this case. 

• (iii) The lightest neutrino is Dirac and has a mass m ~ 3.2 meV. Then, the scalar mass must be 

11.9 meV , ^. 

M = r , the Dirac neutrino contributes. 



If the scalar is the axion, then M = 4.66 meV in this case. 
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Therefore, in all three cases (i)-(iii) where the axion explains the dark energy we predict its mass in the range: 

3.96 meV < M < 4.66 meV . (8.2) 

The left and right ends of the interval in eq. (|8.2[) correspond to no neutrino contribution and to the lightest neutrino 
as a Dirac fermion with mass m — 3.2 meV, respectively. 

In short, we uncovered here the general mechanism producing the dark energy. This mechanism has it grounds on 
well known quantum physics and cosmology. The dark energy appears as a quantum vacuum effect arising when 
stable and weakly coupled quantum fields live in expanding cosmological space-times. No dark energy is produced in 
Minkowski space-time. 

Many research avenues open now connecting dark energy and light particles physics. The more immediate being: (1) 
the study of the radiative corrections to the axion and neutrino cosmological vacuum evolution from their interactions, 
(2) the study of the early neutrino and axion dynamics at temperatures > 1 MeV and > 10 6 GeV, respectively, (3) 
the study of particle propagation in the media formed by the axion and neutrino vacuum, and last but not least (4) 
the probable deep connection between dark energy and dark matter through low energy particle states beyond the 
Standard Model of particle physics. 



APPENDIX A: DIMENSIONAL AND CUTOFF REGULARIZATION OF THE VACUUM ENERGY 



Physical vacuum quantities are computed in sec. IIVI as the equal point limit of two point functions. The distance z 
between the points eci. (|4.1[) naturally plays the role of the regularization parameter. Alternatively, one can regularize 
the two point function with dimensional regularization or cutoff regularization and set z = in the regularized 
expressions. 

In dimensional regularization, we have 

d 4 - 2t k i M 2 ~ 2e 1 



G e (M a, a) =< T tp(x, if) <p{x, rj) >= 



1 



(2 7r) 4 " 2e k 2 -a 2 M 2 + iO (4tt) 2 - £ a 2e 



r(e-i) 



Subtracting the value in Minkowski space-time (a = 1) yields, 



G e (Ma,a)-G e {M,l) = M 2 



[a- 2£ -l] e = u 



M 2 



log a 



(4tt) 2 - £ l j 2(2tt) 2 

in agreement with eq. (|4.8p . 

Alternatively, by regularizing with an ultraviolet cutoff A in four space-time dimensions, we have 



G A (Ma) = \ 



A— >oo 



A 



4ir a 



d 4 k 

(27) 1 P~ 

M 2 



2(2tt) 2 



log 



a 2 M 2 + iO 
A 



Ana Utt/ 



1 + 



A 
Wa 



Ma 



Subtracting the divergence in A = oo again leads to the result eqs. (|A2|) and (|4.8 



G A {Ma) 



A 



4 7r a 



M 2 , 



A 

M 



M 2 



2(2tt) 2 



log a 



(Al) 



(A2) 



(A3) 



We have therefore verified that the point splitting regularization used in sec. IIVI as well as dimensional and cutoff 
regularization methods yield identical results. (It is known since longtime that dimensional regularization gives the 
same physical results than other regularization methods [38|). Analogous results are valid for the two point fermion 
function. 
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